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INTRODUCTION 



The problem of interpretation of scalar mesons is tightly related to tbe most profound topics in particle physics 
which concern the QCD vacuum (see, e.g., the review "Note on scalar mesons" in [l|). These mesons are expected to be 
■ composed of the qq or the lightest 4-quark states, meson-meson molecules or gluonium states. It is disconcerning that 
j up to now a description of this mesonic sector is far from beingcomplete despite of the big effort devoted to studying 
t^. . various aspects of the problem (for recent reviews see, e.g. 0t3). Parameters of the scalar mesons, their nature and 
status of some of them is still not settled [l|. Especially, this concerns the / (600)/cr meson and ifQ(900)/K(800) 
qq , meson. For example, the mass of the former obtained in the Breit-Wigner or iC-matrix approaches ranges in various 
analyses in the interval of about 400-1200 MeV According to the prediction by Weinberg Q based on the mended 
symmetry the mass of the a should be near the mass of the p-meson. As to the mass of the lowest scalar glueball, 
various non-perturbative QCD methods give also very different results. From the QCD sum rules Q one has found 
a scalar-isoscalar meson of the gluonium nature with a mass about 1000 MeV and with the 7T7r-decay width about 
500 MeV. This is in agreement with the recent unquenched-lattice simulation using dynamical fermions Q but it 
diverges from recent calculations on the quenched anisotropic lattices of the glueball spectrum where the mass of the 
lowest glueball is about 1710 MeV [|]. 

The width of the /o(600) (in various experiments and analyses) also has a large spread 600-1000 MeV according 
to an estimate of the Particle Data Group team [l|. Note also the works in which one obtained a very small value 
of 35 ± 12 MeV [l(| and the very large one of about 3200 MeV The prediction for the a- meson width on the 
basis of saturating the superconvergence dispersive sum rules is larger than about 670 MeV [Hj]. The theoretical 
conclusions about widths of glueballs, especially about the lightest one, are also very different in various approaches. 
In Ref. [13j the authors used an effective QCD Lagrangian with the broken scale and chiral symmetry, where a glueball 
is introduced to theory as a dilaton and its existence is related to breaking of scale symmetry in QCD. Then the 7T7t 
decay width of the glueball, estimated using low-energy theorems, is T(G — > tttt) rs 0.6 GeV x (mc / 1 GeV) 5 , where m,Q 
is the glueball mass. I.e., for the glueball with the mass about 1 GeV (if it exists), the width is near 600 MeV. Though 
a use of the above formula is doubtful above 1 GeV, a tendency for the glueball to be wide is apparently seen. This is 
supported by arguments given in [l4| that the glueball width is larger than the ones of the surrounding qq states. On 
the other hand, in Ref. (l5j . where the two-pseudoscalar and two- photon decays of the scalars between 1-2 GeV were 
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analyzed in the framework of a chiral Lagrangian and the glueball was included as a flavor-blind composite mesonic 
held, the glueball was found to be rather narrow in accordance with the former findings of Ref. [l|| . 

Up to now the nature of the /o(980) is not clearly resolved. Besides a qq p~7l. fl8l. [20l f2l | . sub j ect to serious criticism, 
there are recent arguments for a 4-quark state (as the ao(980)) [Hj], a KK molecule [22-24] and a 7777 bound state 
HH23. 

Existence of the /o(1370) meson is still not obvious. In some works, e.g., in [2^, one did not find any evidence 
for the existence of the /o(1370). In Ref. [13] also the best description of irn — > irn, KK was obtained without the 
/o(1370), and it was shown that the KK scattering length is very sensitive to whether this state exists or not. On the 
other hand, in Ref. [3l| a number of data requiring apparently the existence of the /o(1370) is indicated: the Crystal 
Barrel data on pp — > 777771-° [32| and on pp — > 37r° [331, the BES data on JJjh — > </>7r + 7r _ (34|; the / (1370) appears also 
in the GAMS data for 7r + 7r~ — > 7r°7r° at large \t\ [35|. For example, in [36| it was shown within the so-called "hidden 
gauge formalism" that the /o(1370) might be dynamically generated from the pp interaction. 

Especially it is worth to discuss the situation with scalar states in the 1500 MeV region. First, a state, observed 
in this region could be a real candidate for the lightest glueball (see, e.g., Ref. [Hj]). In the mo del- independent 
analyses of data on the processes nn — > nn KK ,r\r\,r\r( using different uniformizing variables [25M27L I30I l37H40l ]. a 
wide state /o(1500) was obtained whereas in many other works, which analyzed mainly the production and decay 
of mesons, as cited in the PDG tables 1]. the rather narrow / (1500) is obtained. Therefore, we have supposed 
[13, Hl| that the wide /o(1500), observed in the multi-channel nir scattering, indeed, is a superposition of two states, 
narrow (qq) and broad (glueball). The former is just observed in the processes of decay and production of mesons. 
An indication about nature of the latter follows from the fact that the /o(1500) is coupled with the approximately 
equal strength with the irir, KK and 7777 systems [HI, H3, [3(3, l37T - l39l | and from the arguments of Ref. [3] on the 



glueball width. These suppositions are in some accordance with the results of the combined K- matrix analysis 41 1 
of the GAMS data on ir~p —> 7r n 7r°n, 777777-, 7777're [42j |. BNL data on n~p —> KKn [Hj] and Crystal Barrel data on 



pp (at rest) 7r 7r 7r ,7r 7r°77,7rV [44j, |45j ], which say that in the 1500 MeV region there are the narrow /o(1500) 
and very wide / (1530±^ ). 

The /o(1710) has most likely the dominant ss component (see, e.g., Refs. [13, EH and the review "Note on scalar 
mesons" in [l[). Note, however, that the QCD sum rules [13] and the K- matrix method (48[ showed that both ,fn(1500) 
and /o(1710) are mixed states with a large admixture of the glueball component. There are also schemes [2g, [49| in 
which the coupling of two gluons (of a scalar glueball) with nn (n is nonstrange u or d quark) appears to be suppressed 
by chiral symmetry [501 ] increasing the relative contribution of the ss component. When assuming this consideration 
to be valid up to energies of the /o(1710), one concludes that this state could be an unmixed glueball [5l| . 

In the scalar-isodoublet sector (except for the well-established state ^(1430)) the possible existence of a very 
broad meson in the 700-950 MeV region is discussed in recent years (see the review "Note on scalar mesons" in [l[). 
E.g., in some recent analyses the authors have found a pole which corresponds to this state Kq (900) [5214581 ] . while no 
such state was seen in the experiment performed by the BaBar Collaboration [59J and in the earlier analyses [60l - [62j . 

In view of all above circumstances, the problems connected with determining the nature of the observed mesonic 
states and their assignment to the quark-model configurations are still open in spite of a large amount of work devoted 
to these problems (see, e.g., Refs. [63H67J). It is clear that resonance parameters should be obtained, if possible, in 
a model-independent way. Here, we present results of the combined three-channel analysis of data on the processes 
7T7T — > 7T7T, KK, 7/7/, 7777' in the channel with the quantum numbers I J PC = + ++ . Study of the Kit scattering in the 
channel with I(J P ) — ^(0 + ) and the role of the strange scalar meson JCq(900) (k(800)) goes beyond the scope of this 
paper and will be discussed in Ref. [68[ . We have used a "model- independent" method [13, H3, H3, [H, 53] based on the 
first principles (analyticity and unitarity) directly applied to the analysis of experimental data. This approach permits 
us to omit a theoretical prejudice in extracting the resonance parameters. It is important that an uniformizing- 
variable method allows to avoid a model dependence when considering resonance contributions. This is possible 
since a main model-independent contribution of the resonance can be given by poles and corresponding zeros on a 
uniformization plane, whereas the possible remaining corrected and mo del- dependent contribution of the resonance is 
supposed to be taken into account in the background. This distinguishes substantially our model-independent method 
from the standard dispersion relation approach based also on analyticity and unitarity where, however, the model 
dependence arises inevitably when saturating dispersive integrals by the contributions of resonances. Then in our 
method, considering the obtained disposition of resonance poles on the Riemann surface, bearing witness to a relative 
strength of coupling with corresponding channels, and resonance masses, we draw conclusions about nature of the 
investigated states. 

Unlike in the previous three-channel analysis of the above processes [25T4271 |39| . in this work we used a new 
uniformizing variable in which we took into account the left-hand branch-point at s = related to the thresholds of 
the 7T7r scattering in crossed channels, in addition to the right-hand branch-points related to the thresholds of the 
analyzed processes. This should diminish considerably dependence of the extracted parameters of resonances on the 
background because the elastic part of the tttt background is stipulated mainly by the contribution of the left-hand 
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cuts. 

The layout of the paper is as follows. In Sec. II we outline the two- and three-coupled channel formalism, where 
pole clusters on the Riemann surface are determined as characteristics of multichannel states and a classification of 
two- and three-channel resonances according to the types of the possible pole-clusters is given. We introduce also the 
new uniformizing variable for the three-channel case taking into account the left-hand branch-point at s — 0, and 
show the disposition of the resonance poles and zeroes related to the various pole-clusters on the unifomization plane 
for the 7T7r-scattering S-matrix element. In Sec. Ill we carry out the combined 3-channel analyses of data on the 
processes tttt — > tttt, KK, r/r/ (variant I) and tttt — ¥ tttt, KK, r\r{ (variant II). In Sec. IV we summarize our conclusions, 
propose an assignment of the scalar mesons lying below 1.9 GeV to lower nonets and discuss the obtained results. 

II. THE COUPLED-CHANNEL FORMALISM IN MODEL-INDEPENDENT APPROACH 

Our model-independent method which utilizes an uniformizing variable can be used only for the 2-channel case 
and under some conditions for the 3-channel one. Only in these cases we obtain a simple symmetric (easily to be 
interpreted) picture of the resonance poles and zeros of the S-matrix on the uniformization plane. The S-matrix 
is determined on the 4- and 8-sheeted Riemann surfaces for the 2- and 3-channel cases, respectively. The matrix 
elements S a p, where a, /3 = 1, 2, 3 denote the channels, have the right-hand cuts along the real axis of the complex s 
plane (s is the invariant total energy squared), starting with the channel thresholds Si (i = 1, 2, 3), and the left-hand 
cuts related to crossed channels. The Riemann-surface sheets are numbered according to the signs of the analytic 
continuations of the quantities s/s — s a as follows: in the 2-channel case 



signs (^ImVs — «i, Im^s - s 2 J = ++, — h , H — (1) 

correspond to sheets I, II, III, IV; in the 3-channel case 

signs ^ImVs - ii, Im\/s - ia, ImVs - S3) =+ + +,- + +, h H h,H , ,- + -, + + - (2) 

correspond to sheets I, II,- • • , VIII, respectively. 

The resonance representations on the Riemann surfaces are obtained using formulas from Ref. [40( (see also Ap- 
pendix), expressing analytic continuations of the S'-matrix elements to unphysical sheets in terms of those on sheet I 
(physical) that have only the resonance zeros (beyond the real axis) around the physical region. Then, starting from 
the resonance zeros on sheet I one can obtain an arrangement of poles and zeros of resonances on the whole Riemann 
surface. 

In the 2-channel case we obtain tree types of resonances described by a pair of conjugate zeros on sheet I only in 
Sn - the type (a), only in S22 ~ (b), and in each of Sn and S22 ~ (c). Then the formulas of the analytic continuations 
of the S'-matrix elements to unphysical sheets [4(| (see also Appendix) immediately give the resonance representation 
by poles and zeros on the 4-sheeted Riemann surface: to the resonances of types (a) and (b), there corresponds a 
pair of complex conjugate poles on sheet III shifted relative to a pair of poles on sheet II and IV, respectively. For 
the states of type (c) one must consider the corresponding two pairs of conjugate poles on sheet III. 

In the 3-channel case we obtain seven types of resonances corresponding to seven possible situations when there are 
resonance zeros on sheet I only in Sn - (a); S22 - (b); S33 - (c); Sn and S22 - (d); S22 and S33 - (e); Sn 
and S33 - (f); Sn, S22, and S33 - (g). Examples for the disposition of poles and zeros, corresponding to some of 
these types of the 3-channel resonances, on the unifomization plane will be given in the next section. 

The resonance of every type is represented by a pair of complex-conjugate clusters (of poles and zeros on the 
Riemann surface). The cluster type is related to the nature of the state. For example, if we consider the tttt, KK 
and 7777 channels, a resonance coupled relatively more strongly to the tttt channel than to the KK and r\r\ ones is 
described by the cluster of type (a). In the opposite case, it is represented by the cluster of type (e) (the state with 
a dominant ss component). The glueball must be represented by the cluster of type (g) (of type (c) in the 2-channel 
consideration) as a necessary condition for the ideal case, if this state lies above the thresholds of the considered 
channels. 

A necessary and sufficient condition for existence of the multi-channel resonance is its representation by one of 
the types of pole clusters. Note that whereas cases (a), (b) and (c) can be simply related to the representation of 
resonances by multi-channel Breit-Wigner forms, cases (d), (e), (f) and (g) are practically lost in the Breit-Wigner 
description. 

One can formulate a model-independent test as a necessary condition to distinguish a bound state of colorless 
particles (e.g., a KK molecule) and a abound state [37], SO, [63j . In the 1-channel case, the existence of the particle 
bound-state corresponds to the presence of a pole on the real axis under the threshold on the physical sheet. In the 
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2-channel case the existence of the bound-state in channel 2 (KK molecule) that, however, can decay into channel 
1 (7T7T decay), would imply the presence of the pair of complex conjugate poles on sheet II under the second-channel 
threshold without the corresponding shifted pair of poles on sheet III. 

In the 3-channel case the bound state in channel 3 (ijij) that can decay into channels 1 (tttt decay) and 2 (KK decay) 
is represented by the pair of complex conjugate poles on sheet II and by the pair of shifted poles on sheet III under 
the r\r\ threshold without the corresponding poles on sheets VI and VII. According to this test, an interpretation of 
the /o(980) as the KK molecule was rejected earlier [Z(|. The reason is that this state is represented by the cluster of 
type (a) in the 2-channel analysis of processes tttt — ¥ tttt, KK and, therefore, does not satisfy the necessary condition 
to be the KK molecule. A further discussion of this topic will be given in the last section. 

Unlike the standard dispersion relation approach in our model-independent method we use an advantage of the fact 
that the amplitude is a one-valued function on the Riemann surface. To this end a uniformizing variable is applied, 
which maps the Riemann surface onto a complex plane. This permit us to use the representation by the pole clusters 
very important for the broad multi-channel resonances. This is impossible in the standard dispersion relation and 
iVmatrix approaches or in the Breit-Wigner one. 

In the combined analysis of coupled processes it is convenient to use the Le Couteur-Newton relations |70j • They ex- 
press the S-matrix elements of all coupled processes in terms of the Jost matrix determinant d(y/ s — s±, • • • , \fs — ijv) 
that is a real analytic function with the only square-root branch-points at ^/ s — s a = 0. The important branch points, 
corresponding to the thresholds of the coupled channels and to the crossing ones, are taken into account in the proper 
uniformizing variable. On the uniformization plane, the pole-cluster representation of the resonance is the good one. 

It is obvious that the main model-independent contribution of resonances is factorized in the S'-matrix elements from 
the background. The possible remaining corrected and model-dependent contributions of resonances are supposed to 
be included in the background. This is realized in a natural way: in the background, the corresponding elastic and 
inelastic phase shifts increase when some channel is opened. Therefore, we denote our approach as model independent. 

In the 2-channel case, the S'-matrix, determined on the 4-sheeted Riemann surface, can be uniformized using, e.g., the 
inverse Zhukovskij transformation [4(| in which the thresholds of two channels are taken into account. Unfortunately, 
already in the 3-channcl consideration a function determined on the 8-sheeted Riemann surface can be uniformized 
only on torus. This is unsatisfactory for our purposes. Therefore, we neglect the influence of the lowest (tttt) threshold 
branch-point at si (however, unitarity on the 7T7r-cut is taken into account). An approximation like this means the 
consideration of the nearest to the physical region semi-sheets of the Riemann surface of the S-matrix. In fact, we 
construct a 4-sheeted model of the initial 8-sheeted Riemann surface approximating it in accordance with our approach 
of a consistent account of the nearest singularities on all the relevant sheets. In practice, neglecting the influence of 
the 7T7r-threshold branch-point means that we do not describe some small region near the threshold. Furthermore, we 
shall take into account in the uniformizing variable also the left-hand branch-point at s — related to the crossed 7T7r 
channels. The allowance for this branch-point should diminish the background dependence of the obtained results. As 
was indicated repeatedly (see, e.g., Ref. [llj) many analyses are subject to criticism (especially with the point of view 
of proof of the resonance existence) because the wide-resonance parameters are strongly controlled by the non-resonant 
background; this is particularly related to low-lying states. The allowance for the left-hand branch-point, related to 
the crossed channels, serves in part for a solution of this problem. For example in Ref. [30( a combined analysis of 
the processes tttt — > tttt, KK in the isoscalar-scalar sector was performed using the method of a uniformizing variable 
which includes two threshold branch-points and the left-hand one at s = 0. In this respect a parameterless description 
of the tttt background was obtained due to the inclusion of the indicated left-hand branch-point. Moreover, it was 
shown that the large background, obtained in earlier analyses of the S-wave tttt scattering [401 ] . hides in reality the 
cr-meson lying below 1 GeV. 

In the 3-channel case the new uniformizing variable used below can have the form: 



w _ \/{s- 52)53 + y/(s- s 3 )s 2 ^ 
a/s(s 3 - s 2 ) 

where we neglect the lowest 7T7r-threshold branch-point and take into account the threshold branch-points related to 
two remaining channels and the left-hand branch-point at s = 0. This variable maps our model of the 8-sheeted 
Riemann surface onto the w-plane divided into two parts by a unit circle centered at the origin. The semi-sheets I 
(III), II (IV), V (VII) and VI (VIII) are mapped onto the exterior (interior) of the unit disk in the 1st, 2nd, 3rd and 4th 
quadrants, respectively. The physical region extends from the point tttt on the imaginary axis (the first tttt threshold 
corresponding to s±) along this axis down to the point i on the unit circle (the second threshold corresponding to 
82)- Then it extends further along the unit circle clockwise in the 1st quadrant to point 1 on the real axis (the 
third threshold corresponding to S3) and then along the real axis to the point b = (y^ + v^3 )/V s 3 — s 2 hito which 
s = 00 is mapped on the w-plane. The intervals (—00, —b], [— b , [b, 00) on the real axis are the images of the 
corresponding edges of the left-hand cut of the 7T7r-scattering amplitude. In Figs[T]and[2J the 3-channel resonances of 
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all the standard types in Sn(w) are represented by the poles (*) and zeroes (o) symmetric to these poles with respect 
to the imaginary axis giving corresponding pole clusters. The "pole-zero" symmetry guarantees the elastic unitarity 
of 7T7T scattering in the (nir, i) interval. For the data we use the results of phase analyses which are given for phase 
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FIG. 1: Uniformization w-plane for the 3-channel-7T7r-scattering amplitude. Representation of resonances of types (a), (b), (c) 
and (d) is shown. 



shifts of the amplitudes d a /3 and for the modules of the S'-matrix elements r] a p — \S a p\ (ot,fl — 1, 2, 3): 

S aa = r] aa e 2lS °o , S af3 = nape**" . (4) 
If below the third threshold there is the 2-channel unitarity then the relations 

f7ll=»722, Vl2 = (1 - »7ll 2 ) 1/2 , 012 = + <?22 (5) 

are fulfilled in this energy region. 
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FIG. 2: Uniformization w-plane for the 3-channel scattering amplitude. Representation of resonances of types (e), (f), and (g) 
is shown. 



Masses and total widths can be calculated using the denominator of formula for the resonance part of the amplitude 
of the form 

T res = 2 ^ rv (6) 

m res - S - ly/sL tot 

taking the pole positions on sheets II, IV and VIII because, as one can see in Ref. [4(j, the analytic continuations of 
the corresponding .S'-matrix elements on these sheets only have the forms 



33 • 



respectively. This means that the pole positions of resonances only on these sheets are at the same points in the 
complex-energy plane as the resonance zeros on the physical sheet and are not shifted due to the coupling of channels. 
Note that the poles on indicated sheets are not always nearest to the physical region. 
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III. ANALYSIS OF THE ISOSCALAR-SCALAR SECTOR 

We analyzed the isoscalar S-waves of the processes 

TTTT — > TTTT, KK, JjJj, rjrj' . 

The experimental data on the tttt scattering from 0.575 GeV to 1.89 GeV are taken from Ref. [7l[ and below 1 GeV 
from the works [72|-|74(. For tttt — > KK the data [75j from threshold to about 1.6 GeV are used. The data for tttt — > rjrj 
up to 1.72 GeV are taken from Ref. [76| and for tttt — > rjrj' from threshold to 1.81 GeV from Ref. [77| . 

In the model- independent approach we have performed two variants of the 3-channel analysis: variant I - the 
combined analysis of processes tttt — > tttt, KK, rjrj; variant II - analysis of tttt — > tttt , K K , 7777' . The influence of the 
^'-channel in the case I and of 7777 in the case II are taken into account via the background. The analysis has been 
carried out with the new uniformizing variable @ (S3 is Am 2 in variant I and (m r] + m,y) 2 in variant II; in the 
following the quantities related to variant II are primed). 

On the w-plane the Le Couteur-Newton relations have the form : 

_ d*(-w*) _ dj-w- 1 ) djw- 1 ) 

, 2 _ d^w*- 1 ) „ „ G2 _ d^-w*- 1 ) Q B G2 _ d{-w) 



S11S22 — S 12 — —, — : , S\\S®& — 5"l3 — "77 — ^ j *5'22S'33 — £"23 — 



d(w) d{w) a{w) 

In this case the subscripts in the matrix elements S a p denote a, f3 =1-tttt, 2-KK, 3 — 7777 or rjrj' . 

The ^-matrix elements in relations ([7]) are taken as the products S = SBS res where Sb describes the background 
and S res the resonance contributions. The d-function for the resonance part is 

M 

d res (w) = w^-^ Y[(w + <) (8) 
r=l 

where M is the number of resonance zeros. For the background part Sb the (i-function has the following form: 

3 / «,_«, 

d B = exp[-i(a + ^ — -(a„ + i0 n ))] (9) 



n=l 



S — So- , S — S 



Sfj Sy 



-O(s-sv), (10) 



s — s s — s. 
/3„ = Kl + Ka -0(s - s a ) + b nv -6(s - s v ) (11) 

S(j Sij 

where s a is the a a threshold, s v the combined threshold of many opened channels in the range of ~ 1.5 GeV 
(17/7', pp, wuS). These threshold are determined in the analysis. 
In variant II, the terms 

, s — Ami , s — Am 2 .-, 

were added to a' n and /3' n , respectively, to account for an influence of the rjrj channel. 

In the analysis, we included all the five resonances discussed below 1.9 GeV in the PDG issue [l|. In variant I, for 
the 7T7r-scattering and tttt — > KK, we considered the data for phase shifts and modules of the S'-matrix elements in the 
energy regions from about 0.4 to 1.89 GeV and from the threshold to about 1.6 GeV, respectively; for tttt — > rjrj, the 
data for the squared module of the S'-matrix element from the threshold to 1.72 GeV. A satisfactory description has 
been obtained. Furthermore it was found 17811 t hat the data on the 7T7T scattering below 1 GeV admit two solutions. 
Let us call them, as in the indicated work [731, " U P" an d "down" because the 7T7r-scattering phase shift goes a bit 
higher in the first case than in the second case. In variant I, for the "up" solution, we considered the representation 
of resonances by different pole-clusters that are admitted by the data. For the "down" solution, we show the formally 
best case. In all cases, the /o(600) is represented by the pole cluster of type (a), the /o(980) is represented only by 
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TABLE I: Variant I: the quality of description of the data for the best variants of representation of considered states obtained 
in the analysis. The letters in the second column denote the pole clusters describing respectively resonances /o(1370), /o(1500) 
and / (1710). 



Solution 




7T7r scattering 
X 2 /dof 




7T7T — )• 

X 2 /dof 




7T7T — > r/TI 

X 2 /ndp 


The total 

X 2 /dof 






bgb 


155.784/(169 - 35) p 


a 1.16 


148.702/(120 - 33) p 


a 1.71 


1.02 


320.770/(305 - 42) p 


a 1.22 


"up" 


cgb 


150.704/(169 - 35) p 


a 1.12 


154.097/(120 - 33) p 


a 1.77 


0.91 


319.447/(305 - 42) p 


a 1.21 




bgc 


146.684/(169 - 35) p 


a 1.09 


152.274/(120 - 33) p 


a 1.75 


0.99 


314.754/(305 - 42) p 


a 1.20 




cgc 


143.269/(169 - 35) p 


a 1.07 


154.469/(120 - 33) p 


a 1.78 


0.91 


312.365/(305 - 42) p 


a 1.19 


"down" 


cgc 


159.425/(169 - 35) p 


a 1.19 


144.475/(120 - 33) p 


a 1.66 


0.73 


315.531/(305 - 42) p 


a 1.20 



the pole on sheet II and shifted pole on sheet III; the resonances /o(1370) and /o(1710) can be described by the pole 
clusters of type (b) or (c); / (1500), of type (g). 

In Table U we demonstrate quality of description for each separate process in the frame of the best combined 
description of all three processes for various acceptable variants of representation of considered states. We use 
abbreviations "dof" - number of degrees of freedom and "ndp" - number of data points. 

It is clear that one should achieve the best description of the separate process; however, then the combined descrip- 
tion of all three processes would be worse. 

In Tables ILTl and UTTl there are given the masses and total widths of states for the indicated cases, calculated from 
the pole positions on sheets II, IV and VIII for resonances of types (a),(b)and (c), respectively, using denominator of 
the resonance part of amplitude in the form (JSJ). For the resonance /o(1500) of type (g), the poles can be used on all 
indicated sheets. 



TABLE II: The "up" solution: the masses and total widths (in MeV) of the /o resonances, obtained at analyzing for acceptable 
variants of representation of considered states. The letters in the upper row denote the pole clusters describing respectively 
resonances /o(1370), /o(1500) and /o(1710). 





bgb 


cgb 


bgc 




cgc 




State 


TYl-res 


Ttot 






TTl-res 


Ltot 


7TL res 


r \ot 


/o(600) 


713.7±5.4 


627.0±7.2 


735.0±6.1 


686.0±7.0 


627.0±7.3 


665.8±11.0 


604.5±5.7 


567.0±5.4 


/o(980) 


1007.6±2.2 


45.2±2.8 


1007.1±2.6 


50.6±2.8 


1007.3±1.9 


50.8±2.8 


1004. 7± 


54.2±2.8 


/o(1370) 


1404.0±7.0 


279.1±22.0 


1390.5±14.3 


223.4±42.8 


1325.6±11.1 


344.6±24.4 


1374.5±16.7 


322.0±60.8 


/o(1500) 


1532.6±15.9 


648.2±26.6 


1544.9±12.2 


646.2±26.0 


1556.6±13.5 


690.4±28.6 


1535.4±12.3 


671.4±26.8 


/o(1710) 


1750.9±35.6 


118.2±30.2 


1751.0±23.8 


118.0±50.8 


1759.2±755.7 


207.0±420.3 


1759.2±716.4 


201.8±385.8 



TABLE III: The "down" solution: the masses and total widths (in MeV) of the fo resonances, obtained at analyzing for the case 
when the resonances /o(1370), /o(1500) and /o(1710) are described by the pole clusters of type (c), (g) and (c), respectively. 



State m res T to t 

/o(600) 769.0±10.0 1036.9±11.8 

/o(980) 1007.2±3.1 64.6±3.0 

/o(1370) 1396.4±24.7 355.2±79.6 

/o(1500) 1534.1±9.2 636.6±25.8 

/o(1710) 1731.0±43.6 203.4±34.8 



As to a quality of description, it is impossible to select any of the above-indicated solutions on the basis of analyzing 
jointly only three considered processes. It is required to add in the combined analysis also relevant processes of 
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decay. We selected the "up" solution mainly because its parameters of the /o(600) remarkably accord with prediction 
(m CT ~ m p and T tot s» 600 MeV) by Weinberg Q, however, for now we should consider both solutions. Furthermore, 
we take a scenario in which the / (1370), /o(1500) and /o(1710) are described by the pole clusters of type (c), (g) 
and (b), respectively. The point is that the parameters of the /o(1500) can be calculated from the pole positions on 
sheets II, IV and VIII. Therefore, an additional criterion for self-consistency of results is a mutual closeness of values 
of the obtained parameters of this important state on the indicated sheets. According to this criterion the selected 
scenario is the most relevant. 

In Tables [TVl and fVl we show the obtained pole clusters for resonances in the complex energy plane y/s, corresponding 
to the cases when the /o(1370), /o(1500) and /o(1710) are described by the pole clusters of the type (c), (g) and (b) 
for the "up" and of type (c), (g) and (c) for "down" solutions, respectively. The poles corresponding to the / (1500) 
on sheets IV, VI, and VIII are of the 2nd order and those on the sheet V of the 3rd order in our approximation. 



TABLE IV: The pole clusters for the /o-resonances for the "up" solution in variant I. y/~Sr~ — E r + iIY/2 in MeV is given. 



Sheet 




II 


III 


IV 


V 


VI 


VII 


VIII 


/o(600) 


Er 

rv/2 


650.1T6.6 
343.0T3.5 


703.8±10.9 
343.0T3.5 






655.9±27.6 
343.0±3.5 


602.2T22.0 
343.0±3.5 




/o(980) 


Er 

T r /2 


1006.8±2.6 
25.3±1.4 


980.8T3.8 
37.2T2.2 












/o(1370) 


E r 

r r /2 








1386.0T14.2 
182.9T34.2 


1386.0T14.2 
179.5±33.6 


1386.0T14.2 
155.7±18 


1386.0T14.2 
111.7±21.4 


/o(1500) 


r r /2 


1510.7±12.2 
323.1±13.0 


1530.0T12.7 
164.7±11.0 


1510.7T12.2 
290.5T46.2 


1510.6T8.5 
156.9T9.0 


1513.T5.8 
193.2±6.4 


1486.4T14.6 
134.4±21.6 


1510.7T12.2 
332.9T73.7 


/o(1710) 


E r 




1750.0T23.8 
65.2±26.6 


1750.0T23.8 
59.0T25.4 


1750.0T23.8 
63.0T23.6 


1750.0T23.8 
69.2T24.0 






TABLE V: The pole clusters for the /o-resonances for the "down" solution in 


variant I. y^iT- = 


E r + T r /2 in 


MeV is given. 


Sheet 




II 


III 


IV 


V 


VI 


VII 


VIII 


/o(600) 


E; 

Er/2 


567.9±12.4 
518.5T5.9 


642.0T17.7 
518.5T5.9 






647.7±29.1 
518.5T5.9 


573.6T25.5 
518.5T5.9 




/o(980) 


E r 

r r /2 


1006.7T3.1 
32.3±1.5 


970.1±5.8 
55.4T2.6 












/o(1370) 


Er 
Tr/2 








1385.1±24.4 
287.0T73.7 


1385.1T24.4 
267.4T83.1 


1385.1T24.4 
158.0T41.8 


1385.1T24.4 
177.6T39.8 


/o(1500) 


Er 
Tr/2 


1500.7±4.9 
318.3±12.9 


1495.0±9.0 
133.6±10.6 


1500.7±4.9 
231.9T17.6 


1496.7T7.2 
141. ±6.3 


1510.2±4.6 
185.2T4.0 


1501.2±9.8 
99.T18.0 


1500.7±4.9 
345.9±14.5 


/o(1710) 


Er 
Tr/2 








1728.0T43.7 
139.9T69.0 


1728.0T43.7 
138.7T8.9 


1728.0T43.7 
100.5T48.4 


1728.0T43.7 
101.7T17.4 



The background parameters are: "up" solution - a = 0.4704± 0.0364, an = -0.2376±0.0132, a Xa = 0.186±0.0335, 
a lv = -0.0788 ± 0.0535, &n = b la = 0, b lv = 0.0305 ± 0.0112, a 2 i = -1.7768 ± 0.0461, a 2a = 0.5204 ± 0.0254, 
a 2v = -9.22±0.649, b 21 = 0.0132 ± 0.0131, b 2a = 0, b 2v = 7.385 ±1.354, b 31 = 0.5494 ±0.0458, b 3a = 0.8995 ±0.0997, 
b 3v = 0; s a = 1.638 GeV 2 , s v = 2.085 GeV 2 ; "down" solution - a = 0.2431 ± 0.0322, a n = -0.0553 ± 0.0113, 
a la = 0.0914 ± 0.0103, a lv = -0.0478 ± 0.0098, b n = b la = 0, b lv = 0.0469 ± 0.0104, a 21 = -1.6811 ± 0.0426, 
a 2a = -0.247±0.1987, a 2v = -7.2±0.5858, b 21 = 0.0329±0.0131, b 2a = 0, b 2v = 7.765±1.4301, b 31 = 0.6135±0.0495, 
6 3ct = 0.6617 ± 0.1099, b 2v = 0. 

The obtained zero positions of the resonances on the Tf-plane are: 
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; 'up" solution - 

for / (600) 

for / (980) 
for / (1370) 

for / (1500) 

wie = 

for / (1710) 



"down" solution 



for / (600) : 

for / (980) : 
for / (1370) : 

for / (1500) : 

w le = 

for / (1710) : 



W! = 3.8538993 + 4.1546473i, w 2 = 

w 3 = -3.8192504 - 4.1137414i, w 4 = 

w 5 = 0.6671266 + 1.1471002i, w 6 = 

w 7 = -3.0340738 - 0.4551088i, w s = 

w 9 = -0.3346148 - 0.0335710i, w 10 

w n = 3.4148479 + 0.5291918z, w 12 

wis = -0.2847651 - 0.0449140i, w u 
w 17 = wig = -3.2842692- 0.3472515i, 

w w = w 20 = 3.2517666 - 0.2894345i, w 21 

w 23 = 0.2818700 + 0.0051870i, w 2i 

W25 = -3.5479990 - 0.0764492i, w 26 



wi = 5.2482964 + 3.3618022i, w 2 = 

w 3 = -4.7510997 - 3.126732LS, w 4 = 

w 5 = 0.7366709 + 1.2022465i, w 6 = 

w 7 = -3.1626888 - 0.6250175i, w 8 = 

w 9 = -0.3228711 - 0.0481947i, w w 

w n = 3.3973448 + 0.5354069i, w 12 

w 13 = -0.2837024 - 0.0468876z, w u 
Wl7 = Wl8 = -3.2730730- 0.336754H, 

wig = w 20 = 3.2 1 6 1 986 - 0.2715142i, w 21 

w 23 = -3.5403847 - 0.1579625i, w 24 

W25 = -0.2829782 - 0.0093685i, w 26 



-0.1320698 + 0.140033H, 

0.1109493- 0.1194408z, 

-0.2290429 + 0.5613030i, 

3.0385254 - 0.4624306i, 
= 0.3351386- 0.0326940*, 
= 0.2899891 + 0.0420199z, 
= Wl5 = -0.3016876 + 0.0264597z, 

= w 22 = 0.31103 6 5 - 0.0259720i, 
= -0.2817805 + 0.0057246i. 
= 3.5470806- 0.0696555i; 



-0.1459568 + 0.0959815i, 

0.1358737 - 0.0873064i, 

-0.2172355 + 0.4780125i, 

3.1950376 -0.6565649z, 
= 0.3267097- 0.0444134i, 
= 0.2981317 +0.0378036i, 
= Wl5 = -0.3096900 + 0.0250295z, 

= w 22 = 0.3114429 - 0.0186598i, 
= 3.5407727- 0.1592757i. 
= 0.2830083- 0.0092620i. 



For variant II we got the following description: for ww scattering x 2 /dof = 148.786/(169 — 30) ~ 1.07; for ww —> KK 
X 2 /dof = 155.006/(120- 29) w 1.70; for ww -> ryn' X 2 / nd P w 0.3. The total X 2 /dof is 306.187/(297- 37) w 1.18. In 
this case the / (600) is described by the cluster of type (a'); / (1370), type (b'); / (1500), type (d'); / (1710), type 
(c'). In Table [VTl wc indicate the obtained pole clusters for resonances on the eight sheets in the complex energy plane 
y/s. The poles on sheets IV and V, corresponding to the /o(1500), are of the 2nd order (this is an approximation). 



TABLE VI: The pole clusters for the /o-resonances in variant II. yfs^, = Ej, + iT' r /2 in MeV is given. 



Sheet 




II 


III 


IV 


V 


VI 


VII 


VIII 


/o(600) 


r;/2 


558.7±13.3 
529T17.4 


564.3T13.7 
529T17.4 






541.3±55.0 
529T17.4 


535.7±26.0 
529T17.4 




/o(980) 


Ej, 

r;/2 


1009.0T3.1 
31.8±1.8 


986.1±5.5 
57.4T2.9 












/o(1370) 


Y' r /2 




1411.6T8.1 
215.6T21.2 


1411.6T8.1 
235T22.6 


1428.4T11.0 
235T22.6 


1428.4T11.0 
215.6T21.2 






/o(1500) 


e', 

r;/2 


1496. 9±4. 7 
198.5±7.8 


1503.0±3.9 
236.0T5.7 


1496. 9±4. 7 
193.1T8.6 


1496.9T4.7 
198.5T7.8 


1494.6±2.9 
193.7T4.7 


1496. 9±4. 7 
193.1T8.6 




/o(1710) 


r;/2 








1743.0T17.8 
144.1T40.3 


1743.0T17.8 
111.5T31.9 


1743.0T17.8 
82.1T36.6 


1743.0T17.8 
114.7T38.6 



The background parameters are: a 1 = 0.2315 ± 0.0085, a' u = 0, a[ v = -0.0616 ± 0.0321, a' la = 0.0298 ± 0.0876, 
a' lv = 0.0622 ± 0.0703, b' u = b[ = b' la = 0, b' lv = 0.0449 ± 0.0105, a' 21 = -3.1359 ± 0.0628, a' 2 = 0, a! 2a = 
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0.4866 ± 0.2778, aL 



&31 



,, - -4.532 ±0.7199, b' 21 = 0, b' 2v 
0.4489 ±0.0606, s a = 1.638 GeV 2 , s v = 2.126 GeV 2 . 
The obtained zero positions of the resonances on the w'-plane are: 



0.7478 ± 0.0607, b' 2a = 2.5545 ± 0.2067, bL = 1-948 ± 1.785 



2r 



for / (600) : 

for / (980) : 
for / (1370) : 

for / (1500) : 

w 15 = 

for / (1710) : 



2.8974465 ± 2.0075214?;, 



II' r-, 



-2.9703646 
0.3793973 4 
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- 2.0251533?;, 
1.1246624?;, 
0.2560839?;, 

- 0.4907520?;, 

- 0.3852687z, 
0.3852687?;, 

= -0.6543233 + 0.1907843?;, w' 17 
w' 19 = -1.5488743 - 0.1115869?;, 
w' 21 = -0.6417003- 0.0474415?;, 



w' 7 = 0.6466646 -) 
w' 9 = -1.3218615 
w' n = 1.3583321- 
= 1.3583321- 



'13 



w 2 



10 



12 



2(1 



= -0.2343026 ±0.1631268?;, 
= 0.2287780- 0.155117H, 
= -0.2267376 ±0.6892463?;, 
= -0.6616945 ±0.2603249*, 
= 1.3490514- 0.5056393i, 
= 0.6855763 ± 0.1928257?;, 
= 0.6855763- 0.1928257z, 
= -1.3502545- 0.3835908i, 
= 1.5602461 - 0.1413653?;, 
= 0.6471581- 0.0350139?;. 



Masses and total widths of the states, calculated from the pole positions on sheets II , IV and VIII for resonances 
of types (a),(b)and (c), respectively, and on sheets II or IV for resonance of type (d), are presented in Table IVTT1 



TABLE VII: Variant II: the masses and total widths of the /o resonances. 



State m' res [McV] T' tot [MeV] 

/o(600) 769.4T15.4 1058.0T34.8 

/o(980) 1009.5T3.1 63.6T3.6 

/o(1370) 1431.0T8.8 469.9T45.2 

/o(1500) 1510.1T4.8 397.1T15.6 

/o(1710) 1746.8T17.9 229.4T77.2 



In Figures [SHSl we show results of fitting to the experimental data in both variants. 

Note, that the "up" solution is not revealed in variant II when in the analysis the data both below and above the 
.ftTiVthreshold are used. This suggests that, the ^-threshold branch-point must be taken into account explicitly. It 
is not sufficient to consider influence of the r\r\ channel only via the background. 

To learn more on existence of the /o(1370) (see the discussion in Introduction), we considered a possibility of 
description in the above-selected cases without this state, i.e. when the resonances /o(600), /o(1500), and /o(1710) 
are represented respectively by the pole clusters of type (a), (g) and (b) in variant I (the "up" solution) and of type 
(a'), (d') and (c') in variant II. The /o(980) is represented by the poles on sheets II and III in both variants. In 
Table [VlTT] we show a quality of description of each separate process for these cases in the frame of the best combined 
description of all three processes. 



TABLE VIII: The quality of description of the data without the / (1370). 



Variant 


■7T7T scattering 


7T7T -> KK 


7T7T — > rjr], r\rj 


The total 




X 2 /dof 


X 2 /dof 


X 2 /ndp 


X 2 /dof 


I 


151.395/(169 - 31) w 1.10 


155.056/(120 - 29) w 1.70 


1.04 


323.14/(305 - 38) ^ 1.21 


II 


150.145/(169 - 26) » 1.05 


160.056/(120 - 26) w 1.70 


0.38 


313.225/(305 - 34) w 1.19 



When calculating x 2 in all cases, the following experimental points have been omitted as obviously strongly falling 
out from the energy dependence: from the 7T7t scattering data the points at 990 MeV for the phase shift Sn and for 
?7ii = 1 5x1 1, and for 77x1 the point at 1650 MeV as strongly violating the unitarity condition. From the 7T7t — > KK 
data there were omitted the points at 1002, 1208.9 and 1235.7 MeV for the 77x2 = |<Si 2 | and the points at 1073, 1082 
and 1387 MeV as giving the anomalously big contribution to x 2 - 
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FIG. 3: The phase shift and modulus of the S'-matrix element in the S-wave 7T7r-scattering. The solid and short-dashed curves 
correspond to variant I, the "up" and "down" solutions, respectively; dash-dotted to variant II. The data are from Refs. [71H74] ]. 



One can see that the description of the 7T7t scattering without the /o(1370) is a bit improved whereas the one of 
the 7T7T — > KK process is made slightly worse, especially as to the phase shift. Generally, an existence of the fo(1370) 
is for now a standard point of view. One ought to take into account also arguments to its favor in Ref. [31| (see 
Introduction). In any case, the existence of the /o(1370) does not contradict the considered data. 

Let us make some more remarks. First, the fact that in variant II we obtain a better description than in variant 
I points to the importance of taking into account the 7777' threshold explicitly. However, in variant II we encounter 
elements of some pseudo-background: these are the negative values of the b coefficients related to an inelastic part of 
the background. The increasing inelastic part of the background implies a necessity to consider explicitly some physical 
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FIG. 4: The phase shift and modulus of the S'-matrix element in S-wave of mr — > KK. The solid and short-dashed curves 
correspond to variant I, the "up" and "down" solutions, respectively; dash-dotted to variant II. The data are from Ref. [75| . 



phenomenon, e.g., additional resonances or representation of resonances by other pole-clusters or the consideration in 
the uniformizing variable of other channel thresholds. The latter situation is the case here: the negative sign of the 
quantity b' 2 = —0.7478 implies the necessity of an explicit consideration of the f^-threshold branch-point. Therefore, 
as to the resonances below 1500 MeV the more adequate description is variant I whereas for the ones above 1500 MeV 
variant II. 

It turns out that the state /o(980) lies slightly above the KK threshold. It is described by the pole on sheet II 
and by the shifted pole on sheet III under the r\r\ threshold without the corresponding poles on sheets VI and VII, as 
it was expected for standard clusters. This may suggest that the /o(980) is not the qq state and can be interpreted, 
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FIG. 5: The squared modules of the irn — >• r/r; (upper figure) and rnr -> r\r( (lower figure) S-wave matrix elements. The data 
are from Ref. [Zg] (upper figure) and from Ref. [73] (lower figure). 



e.g., as a rjrj bound state in accordance with the test discussed in Sec. II: the necessary condition for this is fulfilled. 
See, however, the further discussion of this matter in the last section. 

As to a representation of the /o(600) and /o(980) states, both variants completely agree. The /o(1370) is described 
by the clusters of type (b) or (c) in various scenarios of variant I and of type (b') in variant II; this is reasonable 
taking into account the quark contents of the KK and ryr} systems and the nearness of corresponding thresholds. 
From this we, therefore, deduce that a ss component of the /o(1370) is dominant. This interpretation quite explains 
why one did not find evidence for the existence of the /o(1370) 29] considering only the irir scattering. 

The /o(1500) is described by the cluster of type (g) in variant I and of type (d') in variant II. The former indicates 
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approximately equal coupling constants of this state to the irir, KK and rjrj systems, which apparently could point 
up to its glueball nature. The latter tells on the approximately equal coupling of this state with the irir and KK 
channels whereas the coupling with the rjrj' channel is suppressed; these facts also point up to its glueball nature [l6[ ■ 
Therefore, we deduce a dominant glueball component of the /o(1500). 

Finally, the /o(1710) is described by the clusters of type (b) or (c) in various scenarios of variant I and of type 
(c') in variant II. Taking also into account the quark contents of the rjrj' system this could point to the dominant ss 
component of this state. 

All these conclusions agree quite well with the previous model-independent 2- and 3-channel analyses [25 2 3, H3, H3- 
[40| where other uniformizing variables were used. 



IV. DISCUSSION AND CONCLUSIONS 



The combined analysis of data on the irir — > irir,KK, rjrj, rjrj' processes in the channel with I G J PC = + ++ is 
carried out in the framework of the model-independent approach that is based on analyticity and unitarity and uses 
an uniformization procedure. A new uniformizing variable was used in which, in additional to the right-hand branch- 
points related with the thresholds of the analyzed channels, there is taken into account the left-hand branch-point at 
s = related to the irir scattering in the crossed channels. 

In the analysis of the processes irir — > irir, KK, rjrj it is shown that the data admit two possibilities for parameters 
of the /o(600) with mass, relatively near to the p-meson mass, and with total width about 640 and 1000 MeV. These 
two possibilities are related to two found solutions, admissible by the data below 1 GeV for the phase shift of the irir- 
scattering amplitude: "up" and "down" . As to the combined description of the considered processes it is impossible 
to prefer any of these solutions. However, the "up" solution remarkably accord with prediction by Weinberg [6| with 
respect to the mass and the width, the "down" one to the mass. These values of mass and width, calculated with help 
of formula (HJ) from the pole position on sheet II, correspond to most of the Breit-Wigner values of Refs. [79| (analysis 
of several processes with pseudoscalar mesons) and [80] (GAMS Collaboration, analysis of the reaction pp — > ppir°ir°). 

Furthermore, we have considered all relevant possibilities of representation of resonances by pole clusters (the 3- 
channcl resonances are represented by seven types of the pole clusters). It is shown that for the "up" solution there 
are four scenarios of representation of resonances /o(1370), /o(1500) (as the superposition of two states, broad and 
narrow) and /o(1710) (/o(600) and /o(980) are given by the pole clusters of the same types in all cases) giving about 
the similar description of the above processes and, however, the quite different parameters of some resonances. For the 
/o(600), /o(1370) and /o(1710) a spread of values is obtained for the masses and widths 605-735 and 567-686 MeV, 
1326-1404 and 223-345 MeV, and 1751-1759 and 118-207 MeV, respectively. On the other hancL the results for the 
/o(980) and /o(1500) are more stable and confirm conclusions of our previous analyses [25H27L [30l [37H39I ] . 

Note a quite stable result for the mass and width with rather small errors for the /o(980): m res ss 1005 — 1008 MeV, 
r fot « 45 — 54 MeV. Arrangement of the poles and zeroes on the Riemann surface, which describe this state, may 
suggest that the /o(980) is not the qq state and can be interpreted, e.g., as a rjrj bound state; in any case the necessary 
condition for this is fulfilled. However, following the listings PDG [l(, the mass of this state is obtained above the 
KK threshold in analyses of irir scattering, of multi-channel irir scattering (irir — » irir,KK ,r)r),r]r)') and of processes 
pp(n) —> M1M2M3, whereas below the KK threshold in analyses of the decays of D + -, B + -, J /ip-, and Z-bosons, 
of processes e + e _ — > M1M27, 0M1M27, e + e~ M\Mi, M1M2X, and oi pp —> ppMxM^. Since the mass value below the 
KK threshold is important for a dynamical interpretation of the fo (980) as a KK molecule [22T[24l | it seems that the 
nature of this state is more complicated than a simple rjrj bound state or KK molecule. From the point of view of 
the quark structure these two possibilities are the 4-quark states. It seems this is consistent somehow with arguments 
in favor of the 4-quark nature of /o(980) in work of [19j . 

In view of prolonging discussions of a question, whether the /o(1370) exists or not (see the discussion of this matter 
in Introduction), we considered a description of the multi-channel irir scattering without this state. We concluded 
that an existence of the /o(1370) does not contradict the considered data. The description of the irir scattering is a 
bit improved whereas the one of the irir — > KK process is made worse, especially as to the phase shift. 

The / (1370) (if it exists) and / (1710) have a dominant ss component. Conclusion about the / (1370) agrees 
quite well with the conclusion drawn by the Crystal Barrel Collaboration [sij] where the /o(1370) is identified as rjrj 
resonance in the ir°rjrj final state of the pp annihilation at rest. Interpretation of the /o(1370) as dominated by the 
ss component explains also quite well why one did not find this state considering only the irir scattering. Conclusion 
about the /o(1710) is quite consistent with the experimental facts that this state is observed in 77 — > KsKs [HJ but 
not observed in 77 — > ir + ir~ [83| . 

As to the /o(1500) (m re s = 1510 MeV, r tot = 397 MeV) we suppose that it is the eighth component of octet mixed 
with the glueball being dominant in this state. Its largest width among the enclosing states points also to its glueball 
nature [8J|. Note that in the PDG tables on the /o(1500) listing, an average value for the width of 109 ± 7 MeV 
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is cited. However, there one indicates only the results of analyses of meson production processes, and in the few 
cases where the results of combined analyses of coupled processes are cited, authors did not use the representations of 
the multi-channel resonances by pole clusters (this is especially important in the case of wide resonances), i.e., they 
did not apply all aspects of the multi-channel analysis. On the other hand, one can see from the data on scattering 
processes, analyzed here [7l| . that the energy dependence of observed quantities do not demonstrate a pronounced 
structure in the 1500 MeV region, which is needed for the narrow resonance. Therefore, it is reasonable to suggest 
that in this region there is a superposition of two states, a wide and a narrow one. 

It is known that there is a number of properties of the scalar mesons which do not allow for a satisfactory setup the 
lowest nonet. The main observations are the approximate equal masses of Jo (980) and oo(980) and the ss dominance 
in the wave function of the /o(980). If these states are in the same nonet then the /o(980) must be heavier than 
a (980) by about 250-300 MeV due to the mass difference of s- and u-quark. Exclusion of the / (980) as a non qq 
state and discovery of the -Kfg-doublet (if it will be confirmed) moves off a number of these problems. 

One can propose the following assignment of scalar mesons lying below 1.9 GeV to lower nonets [27|. The lowest 
nonet: the isovector ao(980), the isodoublet Kq(900), and /o(600) and /o(1370) as mixtures of the 8th component of 
octet and the SU(3) singlet. Then the Gell-Mann-Okubo (GMO) formula of 

3m 2 /8 = K, - m 2 aa (12) 

gives m/ 8 = 910 MeV. For this nonet is seems to be important to test the nature of strange scalar meson JCq(900) in 
a model-independent way. This will be the subject of a forthcoming paper [68| . 
In the relation for the masses of the nonet 



m a + w/ (i370) = 2mx« (13) 

the left-hand side is by about 18 % larger than the right-hand one. 

For the next nonet of the radial excitations we find: a (1450), ^(1450), and / (1500) and / (1710), the / (1500) 
being mixed with a glueball which is dominant in this state. From the GMO formula we set to/ 8 w 1453 MeV. In 
the formula 



W/o(1500) + ^/ (1710) = 2mj f *( 14 50) (14) 

the left-hand side is by about 12.5 % larger than the right-hand one. 

This assignment removes a number of prior questions and does not rise new ones. The mass formulas indicate to a 
non-trivial mixing scheme. Breaking of the relations (fT"3|) and (|14[) tells us that the a— /o(1370) and /o(1500)— /o(1710) 
systems get additional contributions absent in the i^g(900) and K g (1450), respectively. A search of the adequate 
mixing scheme is complicated by the circumstance that here there is also a remainder of chiral symmetry, though, on 
the other hand, this permits one to predict correctly, e.g., the cr-meson mass. 
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Appendix A: Analytic continuation of the 3-channel S-matrix elements to unphysical sheets 

Here we show, for convenience, formulas of the analytic continuations of the 3-channel S-matrix elements to 
unphysical sheets of the Riemann surface in terms of those on sheet I (the physical sheet) - S^r that have only 
zeros (beyond the real axis) corresponding to resonances, at least, around the physical region. In Ref. [4(| the general 
formula was given for the case of N channels and as example for three channels. The direct derivation of these formulas 
requires rather bulky algebra. It can be simplified if we use a circumstance that the i^-matrix has the same value in 
all sheet of the Riemann surface of the S-matrix. This fact follows from Hermiticity of the i^-matrix K = K + which 
means that the iC-matrix does not try discontinuity when going across unitarity cuts. Then after some algebra, one 
can obtain formulas under interest shown below in the table. In Table IIX1 the superscript I is omitted to simplify 
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TABLE IX: Analytic continuations of the 3-channel S-matrix elements to unphysical sheets 



Process 


I 


II 


III 


IV 


V 


VI 


VII 


VIII 


1 -> 1 


Sn 


I/S11 


S22/-D33 


D 33 /S22 


det S/Du 


Dn/detS 


S33/D22 


-D22/S33 


1 -> 2 


Sia 


iSi2/Sn 


— S12/-D33 


IS12/S22 


iD 12 /D n 


-D 12 /detS 


iDi 2 /D 2 2 


D12/S33 


2^2 


S22 


D 33 /Sn 


Su/Dsa 


1/S22 


S 33 /Dn 


D 22 /detS 


det S/D22 


D n /S 33 


1 -> 3 


S13 


iSis/Sn 


-iD 13 /D 33 


—D13 1 S22 


-iDi 3 /Du 


Dn/detS 


-S13/D22 


iSl 3 /S33 


2^3 


S23 


Daa/Su 


iD 23 /D 33 


iS2 3 / ' S22 




-D 23 /detS 


iD 2 3/D 2 2 


1S23/S33 


3^3 


S33 


D 22 /Sn 


det S/D 33 


Du/Ska 




D 33 /detS 


Sn/D 2 2 


I/S33 



the notation, det S is the determinant of the 3x3 ^-matrix on sheet I, D a p is the minor of the element S a p, that is, 

Dn — S22S33 — 2 3, D22 — S11S33 — S 13 , D33 = S11S22 — S 12 , D\2 = S12S33 — S\3D23, ^23 = S11S23 ~ S12S13, etc. 
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